Solitary waves on vortex lines in Ginzburg— Landau models 
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Axisymmetric disturbances that preserve their form as they move along the vortex lines in uniform 
Bose- Einstein condensates are obtained numerically by the solution of the Gross-Pitaevskii equation. 
A continuous family of such solitary waves is shown in the momentum (p) - substitution energy 
(£) plane with p — > 0.09pft 3 /c 2 , £ — > 0.091pK 3 /c as U — > c, where p is the density, c is the speed of 
sound, k is the quantum of circulation and U is the solitary wave velocity. It is shown that collapse 
of a bubble captured by a vortex line leads to the generation of such solitary waves in condensates. 
The various stages of collapse are elucidated. In particularly, it is shown that during collapse the 
vortex core becomes significantly compressed and after collapse two solitary wave trains moving in 
opposite directions are formed on the vortex line. 
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An important role in the dynamics of nonlinear sys- 
tems is played by solitary waves - the localised dis- 
turbances of the uniform field that are form-preserving 
and move with a constant velocity. They appear in di- 
verse contexts of science and engineering, such as fluid 
dynamics, transport along macromolccules, fibre optic 
communications just to name few. Considerable inter- 
est attaches to determining the entire sequence of soli- 
tary waves as they define possible states that can be 
excited in the system. Understanding the production, 
motion and interactions of such solitary waves is one of 
the most significant questions in nonlinear science. In 
condensed matter systems solitary waves are topological 
objects since they owe their existence and perseverance 
to the topology of the order parameter field describing a 
medium with a broken symmetry. In this Letter I estab- 
lish and study the production of a new class of solitary 
waves that move along vortex lines/topological defects in 
conservative Ginzburg-Landau systems. The discussion 
will be restricted to condensed matter system such as 
atomic Bose-Einstein condensates (BEC) where the evo- 
lution equation is the Gross - Pitaevskii (GP) model 0; 
see equation below (also known as defocusing non- 
linear Schrodinger equation in nonlinear optics). The 
applications are not restricted to the condensed matter 
systems due to the generality of the Ginzburg-Landau 
systems with implications to the motion of excitations 
along cosmic strings in the early Universe [2| and along 
topological defects in other ordered media: liquid crystals 
Q, non-equilibrium patterns, etc. Finally, the new so- 
lutions to the nonlinear Schrodinger equation which this 
Letter presents are of interest when so few have been 
derived in multidimensions. 

Symmetry-breaking transitions in equilibrium systems 
can be described by an energy functional, £ — J C dV, 
in terms of a Lagrangian L. The simplest form of the 
Lagrangian, capable of describing the quenching beyond 
the critical point where the disordered state becomes 
unstable and the symmetry is spontaneously broken, is 
a Ginzburg-Landau Lagrangian C — ^[IVi/"! 2 + ^(1 — 



|^| 2 ) 2 ] dependent on a complex scalar order parameter 
The evolution equation describing the relaxation to 
the equilibrium state in an energy-preserving (conser- 
vative) system is given by the Euler-Lagrange equation 
tpt = —\dE/dip* that we write as 



v 2 ip + (i - M 2 )V>. 



(i) 



A phase singularity of a complex field ip given by ip = is 
called a quantised vortex or topological defect depending 
on a particular application. The total change of phase 
around any closed contour must be a multiple of 2tt and 
only quantised vortices with the total change of phase 2tt 
are topologically stable. 

Such a general reasoning gives a simple explanation 
why equation Q has a universal meaning and has been 
applied to a variety of systems Q. In particular, it de- 
scribed accurately both equilibrium and dynamical prop- 
erties of BEC The GP model has been remarkably 
successful in predicting the condensate shape in an exter- 
nal potential, the dynamics of the expanding condensate 
cloud and the motion of quantised vortices; it is also a 
popular qualitative model of superfluid helium. For these 
systems equation |QJ is written in dimensionless variables 
such that the unit of length corresponds to the healing 
length £, the speed of sound c = l/y2, and the density 
at infinity = |i/>oo| 2 = L 

The straight-line vortex positioned along the z— axis 
is obtained by rewriting Q in cylindrical coordinates 
(s, 8, z) and using the ansatz ipo = R(s) exp(i0). The 
resulting solution for R(s) was found numerically in the 
first reference of 0] and approximated in 0. The in- 
finitesimal perturbations of a rectilinear vortex in the GP 
model may be bound or free, depending on their angular 
and axial wavenumbers, m and k. The free waves radiate 
energy acoustically to infinity, while the bound states do 
not. The low frequency modes m = 1, that displace the 
axis of the vortex, are found to be bound for all k 0, • 
The low frequency modes m = 2 are also bound when k 
is sufficiently large, but are free for small k 0. 
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A family of fully three-dimensional solitary waves was 
found by Jones and Roberts [j| who integrated nu- 
merically and determined the entire sequence of soli- 
tary wave solutions of the GP equation, such as vortex 
rings and finite amplitude sound waves named rarefac- 
tion pulses. They showed the location of the sequence 
on the momentum, p, energy, £ , plane, that I shall re- 
fer to as the JR dispersion curve. In three dimensions 
they found two branches meeting at a cusp where p and 
£ assume their minimum values, p m and £ m . As p — ► oo 
on each branch, £ — > oo. On the lower branch the solu- 
tions are asymptotic to large vortex rings. As £ and p 
decrease from infinity along the lower branch, the solu- 
tions begin to lose their similarity to large vortex rings. 
Eventually, for a momentum po slightly greater than p m , 
they lose their vorticity (ip loses its zero), and thereafter 
the solitary solutions may better be described as 'rar- 
efaction waves'. The upper branch consists entirely of 
these and, as p — > oo on this branch, the solutions ap- 
proach asymptotically the rational soliton solution of the 
Kadomtsev-Petviashvili Type I equation. 

In what follows I determine the entire family of ax- 
isymmetric (m = 0) solitary wave solutions that move 
along the straight-line vortex and relate them to the JR 
dispersion curve. In equation Q written in cylindri- 
cal coordinates (s, 8, z) I take the ansatz ip = (R(s) + 
<j)(s, z)) exp(i6), and assume that the disturbance <f>(s, z) 
moves with velocity U in the positive z— direction. In the 
frame of reference moving with the solitary wave, <f>(s, z) 
satisfies 



and the energy of the vortex line, 
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The disturbance is localised, so the boundary condition 
is </>(s, z) — > 0, as |x| — > oo in all directions of x. In view 
of the asymptotic expansions at infinity U , I introduce 
stretched variables z' = z and s' = syl — 2U 2 and map 
the infinite domain onto the box (0, ^) x (—■§,§) using 
the transformation z = tan -1 (Lz') and s = tan _1 (Ls'), 
where L is a constant ~ 0.1 — 0.4. Transformed equa- 
tion was expressed in second-order finite difference 
form using 250x200 grid points, and the resulting non- 
linear equations were solved by Newton-Raphson itera- 
tion procedure using banded matrix linear solver based 
on bi-conjugate gradient stabilised iterative method with 
preconditioning. For each solitary wave two quantities 
were calculated: the nonzero (z—) component of the mo- 
mentum j^] 
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and the substitution energy, £, which is the difference 
between the energy of the vortex-solitary wave complex 



£ = \\ |v^| 2 -|VV>ol 2 + ^(i-H 2 ) 2 4(H<Aol 



'dV. (4) 



By performing the variation xjj — ► ^ + 8ij) in © (01 and 
discarding surface integrals that vanish provided Sip — > 
for x — > oo, we see that U = d£/dp, where the deriva- 
tive is taken along the solitary wave sequence. The same 
expression is obeyed by the sequences of classical vortex 
rings in an incompressible fluid and by the solitary waves 
of |9j. I also note that if we multiply £ by zdifj*/dz and 
integrate by parts we get £ = J \dip/dz\ 2 dV, similarly to 
the expression for the energy of the JR solitons 10]. This 
expression was used as a check of the numerical accuracy. 
As p — > oo, £ — ► oo and solitary wave solutions are repre- 
sented by large vortex rings moving along the vortex line. 
As £ and p decrease from infinity, the radii of the rings 
decrease and for a momentum po < 78, </>(s, z) = on 
z— axis only. To distinguish these solutions from vortex 
rings and to emphasise the analogy with the JR solitary 
waves, these solutions will be called rarefaction waves as 
well. Table 1 shows the velocity, substitution energy, mo- 
mentum and radius of the solitary wave solutions found. 
Figure 1 shows the momentum-energy curve of the solu- 
tions in comparison with the JR dispersion curve. Notice 
that unlike the JR dispersion curve, there is no cusp on 
the energy-momentum plane. As U — > c neither £ nor p 
go to infinity, instead £ —> 32 and p — > 45 which lies below 
the JR cusp. The substitution energy, £, and momentum, 
of our vortex rings are larger that corresponding values 
of £ and p of the JR rings moving with the same velocity. 
If the vortex rings of the same radii are compared, our 
rings have lower energy and momentum. Figure 2 shows 
the density isoplots, the density contour plots, and the 
velocity vector fields of two representative solutions such 
as a rarefaction pulse and a vortex ring. Notice the vor- 
tex core expansion at the centre of the vortex ring due 
to a decrease in pressure in high velocity region. 

Table 1. The velocity, U, substitution energy, £, momentum, 
p and radius, b, of the solitary wave solutions moving along the 
straight-line vortex. 



u 



0.4 0.45 0.5 0.55 0.6 0.61 



142 113 90.7 72.4 56.9 54.0 
262 193 145 110 83.2 78.4 
4.18 3.62 3.08 2.41 1.05 0.1 



U 



0.63 0.65 0.67 0.69 0.7 0.75 



48.4 43.0 37.8 33.3 32.2 32.1 
69.4 61.0 53.1 46.5 45.0 45.0 



A question that arises after new solutions are found 
theoretically is how to create them in an actual physical 
system. In jjjj we established a new mechanism of vor- 
tex nucleation by collapsing bubbles in the context of the 
GP model. These results referred to collapse of cavitated 
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FIG. 1: (color online) The dispersion curves for two families 
of the axisymmetric solitary wave solutions. The dashed line 
represents the JR dispersion curve. The part of the curve 
that corresponds to the vortex rings is shown in grey (red). 
The solid line gives the substitution energy as a function of 
momentum for the solitary waves moving along the vortex 
line with vortex rings shown in light grey (green). 
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bubbles generated by ultrasound in the megahertz fre- 
quency range that have been observed to produce quan- 
tised vortices in superfluid helium [l^. Also, vortices 
form as a result of bubbles colliding during a first-order 
phase transition of an early Universe we have 

shown that a soft bubble, carved out in the surround- 
ing fluid by an electron through its zero-point motion, 
becomes trapped in vortex lines. The Bernoulli effect of 
the flow created by the flow circulation around the vortex 
propels the bubble and vortex towards one another with 
a force approximately proportional to s~ 3 , where s is the 
closest distance between them. As the bubble becomes 
trapped in the vortex core, the flow round the bubble 
acquires circulation that it previously could not posses. 
After the emission of Kelvin waves, that were excited 
on the vortex core during the capture, the bubble- vortex 
complex stabilises to an axisymmetric form depicted in 
Figure 2 of see a ^ so Figure 3 (t = 0) below. One 
could expect that a similar capture of bubbles created 
by ultrasound takes place in BEC and that vacuum bub- 
bles get trapped in cosmic strings. The captured bubble 
will than collapse sending axisymmetric waves along the 
vortex line. To elucidate the stages of this collapse I per- 
formed the numerical simulations of the GP equation JJJ 
starting with the initial condition 



V(x,t = 0) 



R(s)e ie ttmh(^) 




if r > a 

if < r < a, 



(5) 



where r 2 = x 2 +y 2 + z 2 . Initial state (JSJ gives an accurate 
representation of the stationary complex, that consists of 
the straight-line vortex at x = 0, y = and the bubble of 
radius a centred at the origin. The surface of the bub- 
ble is assumed to be an infinite potential barrier to the 
condensate particles, so no bosons can be found inside 



FIG. 2: (colour online) Two solitary wave solutions moving 
with velocity U = 0.45 (left side) and U = 0.69 (right side) 
along the straight-line vortex. The density isoplots at p = 0.1 
(left) and p = 0.3 (right) are shown on two top panels. The 
density contour plots (0.1, 0.3, 0.5, 0.7, 0.9) at the cross-section 
9 — are shown on two middle panels. The velocity fields with 
contour plots (0.05,0.1,0.3,0.5) of the velocity magnitude, 
\/v'i + v'i , are given at the bottom. 

U = 0.45 U = 0.69 
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the bubble (-0 = 0) before the collapse, and this is why 
tanh(£/v / 2), which is the wavefunction of the conden- 
sate distance £ away from a solid wall, is relevant here. 
'Softer' bubbles that allow some condensate penetration 
were also considered by reducing the slope of the hyper- 
bolic tangent, but no significant difference was detected. 
I performed fully three-dimensional calculations for cav- 
ities of various radii in a computational cube with sides 
of 200 healing lengths [l4|. Figure 3 shows the density 
plots of the portion of the cross-section at y = at vari- 
ous times after the collapse of a bubble of radius a = 10. 
Figure 4 depicts density contours p/poo = 4/5 at times 
t = 0, 20 and 60. The time-dependent evolution of the 
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FIG. 3: (colour online) Collapse of the bubble of radius a = 10 
trapped by the vortex line. The density plots of the cross- 
section of the solution of Q with initial state © at y — 0, x € 
[—50,50], z G [0,100] are shown. The insert shows a vortex 
ring travelling along the z— axis generated after collapse of 
the bubble of radius a — 50 at t = 100. Two vortex rings of 
smaller radii are also seen on the left and on the right of it. 
Both low and high density regions are shown in darker shades 
to emphasise intermediate density regions. 
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condensate during and after bubble's collapse involves 
several stages. During the first stage dispersive and non- 
linear wave-trains are generated at the surface of the bub- 
ble. This stage of the evolution is characterised by a flux 
of particles towards the centre of the cavity. This creates 
an inward force acting on the vortex core reducing the 
cross-sectional area of the core. The reduction by a fac- 
tor of 1.5 is seen on t = 15 snapshot of Figure 3 and on 
the density contour of Figure 4 at t = 20. The next stage 
in the evolution is outward expansion of the condensate 
that overfilled the cavity. The instability mechanism for 
collapsing bubbles in the absence of the straight-line vor- 
tex, which we described in detail in 0], sets in, leading 
to the production of vortex rings and rarefaction pulses 
mostly along the vortex line, as the energy and momen- 
tum necessary for their creation is lower there. As the 
train of solitary waves starts moving away from the col- 
lapsed bubble (see Figures 3 and 4), the distance between 
them increases since they move with different velocities. 
In time each individual solitary wave approaches its lo- 
calised form found in the first part of this paper; see the 
insert of Figure 4. During collapse of a bubble of a larger 
radius (a > 28) vortex rings will be generated together 
with rarefaction pulses on the vortex line; see the insert 
of Figure 3. 

The author is grateful to Professors Paul Roberts 



FIG. 4: (colour online) Collapse of the bubble of radius 
a = 10. The contour p/ Poo = 4/5 is shown at t = (green 
solid line), t = 20 (blue dotted line), and t = 60 (red dashed 
line). The density contour plot of the condensate showing a 
well-separated solitary wave at t — 166 is given in the insert 
(compare this with the middle right panel of Figure 2). 
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